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The moon-forming impact and the subsequent evolution of the proto-Earth is strongly dependent 
on the properties of materials at the extreme conditions generated by this violent collision. We 
examine the high pressure behavior of MgO, one of the dominant constituents in the earth’s mantle, 
using high-precision, plate impact shock compression experiments performed on Sandia National 
Laboratories Z-Machine and extensive quantum simulations using Density Functional Theory (DFT) 
and quantum Monte Carlo (QMC). The combined data span from ambient conditions to 1.2 TPa and 
42,000 K, showing solid-solid and solid-liquid phase boundaries. Furthermore our results indicate 
under impact that the solid and liquid phases coexist for more than 100 GPa, pushing complete 
melting to pressures in excess of 600 GPa. The high pressure required for complete shock melting 
places a lower bound on the relative velocities required for the moon forming impact. 

PACS numbers: 91.60.Hg, 91.60.-x, 81.40.Vw, 71.15.Pd 


The leading theory of moon formation is a giant im¬ 
pact event occurring approximately 4.5 x 10 9 years ago 
[l|. Recent large scale simulations show that giant im¬ 
pacts are capable of producing the correct Earth - Moon 
masses and angular momentum 0I- Complicating the 
giant impact theory, however, is that the Earth and moon 
have a nearly identical chemical and isotopic composi¬ 
tion 0. This implies that either the impactor was com- 
positionally similar to the proto-Earth, which is thought 
to be unlikely, or that extensive mixing of the post impact 
materials occurred. The simulations needed to test the 
post-impact mixing are dependent on an accurate under¬ 
standing of the properties of mantle materials at extreme 
pressures and temperatures. For example, post impact 
mixing for chemical equilibration in the proto-Lunar disk 
has been shown in simulations 0, but requires melting 
and vaporization of the mantle in order for material to 
diffuse. Testing such theories is difficult because the melt 
line of the most common mantle materials is not well con¬ 
strained at these conditions 0. 


Advanced facilities for performing dynamic materials 
experiments have greatly increased the pressure and tem¬ 
perature regimes that can be probed for important plane¬ 
tary materials EM- The Rankine-Hugoniot (RH) equa¬ 
tions [11] relate the experimentally measured shock veloc¬ 
ity to the thermodynamic state, pressure and density, of 
the shocked material. The ability to perform experiments 
with steady planar shocks and with well-characterized 
impactors and targets is critical for determining the equa¬ 
tion of state (EOS) and the phase. To fully address 
the physics relevant to planetary science, this thermo¬ 
dynamic information must be augmented with an under¬ 
standing of the phase transformations. 


In this work we focus on MgO, the end member of mag- 
nesiowristite, a major constituent of the earth’s mantle 
0 and likely other terrestrial exoplanets 0. At ambi¬ 
ent conditions, MgO exists in a NaCl (Bl) lattice struc¬ 


ture, which is stable over a wide pressure-temperature 
13 - 0 . mic compression experiments start- 


range 

ing from ambient Tt]- 21| and from To=1900 I< 221 show 
no indications of phase transitions up to 230 GPa. Ab ini¬ 
tio studies on the MgO phase diagram agree that three 
phases exist: the Bl solid, the B2 (CsCl) solid, and the 
liquid 230, but disagree on the location of the phase 
boundaries. Along the Hugoniot, which is relevant for 
planetary impact scenarios, the location of the B1-B2 
transition and the melt transition has not been precisely 
determined. Recently McWilliams et al. have shown 
that MgO can be dynamically compressed to pressures 
exceeding 1 TPa using a decaying shock technique [ 5 ( 1 ] . 
The authors also proposed locations for the B1-B2 and 
B2-liquid transitions along the Hugoniot, but the mea¬ 
surements relied heavily on prior MgO Hugoniot data, 
which was not well known above 200 GPa. 


We present a comprehensive study of the MgO Hugo¬ 
niot using experiments, density functional theory (DFT), 
and quantum Monte Carlo (QMC) methods over a wide 
pressure range covering the Bl, the B2, and the liquid 
phases from 0.27 to 1.2 TPa. The high-precision data 
constrain the Hugoniot at multi-Mbar pressures, and the 
DFT and QMC results further elucidate information on 
the phase boundaries. This work provides accurate EOS 
data at extreme conditions and furthermore reveals lower 
limits of the relative velocity required in the giant impact 
scenario for moon formation. 

To attain planetary impact conditions, we performed a 
series of shock compression experiments using the Sandia 
Z-Machine [ 53 ] ■ The Z-machine is a pulsed power system 
capable of producing shaped current pulses and induced 
magnetic fields in excess of 20 MA and 10 MG respec¬ 
tively. The combined current and magnetic field densities 
generate magnetic pressures up to 650 GPa that can ac¬ 
celerate aluminum flyers up to 40 km/s [ 53 ]. 

Figure [l] shows a schematic view of the target geom- 
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FIG. 1. The experimental target and representative VISAR 
data. The VISAR measures the flyer velocity as it approaches 
the MgO (grey line). For this low velocity impact, as the 
shock transits the MgO, the VISAR loses signal (red line). As 
the shock transits into the quartz, the VISAR signal returns 
and the quartz shock velocity is measured (blue line). The 
inset shows a VISAR raw signal. The fiducials correspond to 
impact and to shock transit into the rear window. 


etry; a more detailed Z target geometry is found else¬ 
where [29[. An A1 or Cu flyer plate is shocklessly acceler¬ 
ated toward the target stack consisting of a single-crystal 
MgO sample ([100], 300-500 pn i, Asphera Corp., p 0 = 
3.584 g/cm 3 ) backed by a window. Although the back 
side of the flyer is melted by the high current, the surface 
of the flyer impacting the sample typically has a few hun¬ 
dred microns of solid density material at impact [28j. The 
velocity interferometer system for any reflector (VISAR) 
measures the flyer plate velocity (Vf) up to impact at 
the target (Fig. [TJ. Impact produces a steady shock in 
the MgO sample. At low impact velocities (and conse¬ 
quently low shock pressures), the MgO sample scatters 
light from the VISAR preventing direct measurement of 
the shock velocity. Instead, sharp fiducials are observed 
in the VISAR signal (Fig.[T]inset) that correspond to im¬ 
pact and to shock transit into the quartz window. In this 
case, we calculated the MgO shock velocity (Us) using 
the transit time determined from the fiducials and the 
measured thickness. At high impact velocities (and high 
shock pressures), the shock front in the MgO is reflective 
and the VISAR directly measures the shock velocity. As 
the shock transits into the quartz, the VISAR measures 
the quartz shock velocity directly. Multiple VISAR sig¬ 
nals were recorded for each sample eliminating 27 t am¬ 
biguities and providing redundant measurements for im¬ 
proved precision. For directly measured velocities, the 
uncertainty is better than 1% and for transit time mea¬ 
surements the uncertainty is on the order of 1-2%. 

Knowing the initial densities of the MgO and the flyer 
plate and measuring the flyer velocity and the MgO 



FIG. 2. The MgO Hugoniot data in p — P space from the Z 
experiments, previous experimental data [l7H2ll . [26| , and the 
DFT simulation results from this work. The current data set 
shows a deviation from the extrapolation of the fit to the B1 
MgO data from < 230 GPa. 



FIG. 3. The experimental Us —Up data including results from 
previous experiments from Refs [TtI [2~il |26| . The optimized 
linear fits are also plotted. The dashed vertical lines indicate 
the optimized phase boundaries along the principal Hugoniot 
and the shaded cyan regions indicate the uncertainty. 

shock velocity, we calculate the MgO Hugoniot state 
density (p), pressure (P), and particle velocity (Up). 
The Hugoniot state is determined using a Monte Carlo 
impedance matching analysis M to solve the RH equa¬ 
tions a The Monte Carlo method accounts for the 
uncertainties in the experimental measurement and the 
A1 and Cu Hugoniot standards. Further discussion on 
the impedance matching method along with the list of 
Hugoniot data are found in the supplement j30| . 

Figure [2] plots the experimental and DFT principal 
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Hugoniot in p -P space. The Z experimental data span 
the range from 0.27 TPa up to 1.2 TPa - the highest, di¬ 
rectly measured Hugoniot states attained in MgO. Also 
included are the DFT simulation results for the Bl, B2, 
and liquid phases of MgO (discussed later). Although 
the VISAR diagnostic does not give direct information 
about the MgO phase upon shock compression, we can 
infer phase transitions given our data. Figure [2] shows 
an extrapolation of the linear fit to the Us — Up data for 
B 1-phase Hugoniot states < 230 GPa (converted to p- 
P using the RH equations) determined from the previous 
experiments [1711211. Below rs 360 GPa, the Z experimen¬ 
tal data are consistent with the gas-gun data but above 
360 GPa they deviate from the extrapolation. This sug¬ 
gests that the Bl phase is stable up to 360 GPa and likely 
undergoes a phase transition from the Bl state to another 
phase, presumably the B2 state, at that shock pressure. 
At higher pressures («>700 GPa) the shock front in the 
MgO was reflective. From the reflective shock front, we 
infer that the MgO has melted into a conductive fluid, 
similar to what is observed for quartz [7] . 

The interpretation for the MgO phase along the Hugo¬ 
niot is not as obvious between the low pressure Bl phase 
and the high pressure liquid phase. Prior ab initio calcu¬ 
lations suggest that MgO transitions from Bl to B2 and 
then from B2 to liquid with a small B1-B2 coexistence 
region and a potentially large coexistence region between 
the B2 and liquid phases |25[. To investigate the phase 
regions, we analyze the Us—Up data using a Monte Carlo 
optimization (MCO) method similar to the method used 
in previous work on carbon [3 lj. 

Slope changes in the Us — Up data often indicate 
phase transitions and phase boundaries. Using the MCO 
method, we fit four lines to the experimental Us — Up 
data. The Us — Up data were converted to a “cloud” of 
points (as described in Ref. EH), allowing region bound¬ 
aries to move smoothly during optimization. For a par¬ 
ticular set of data clouds, the eleven parameters (four 
slopes, four intercepts, and three region boundaries) of 
a four-segment curve were obtained by minimizing the 
square minimum distances to each cloud point. Revised 
clouds were generated by randomly drawing a new cen¬ 
ter for each cloud. Optimization was repeated (« 10000 
times) using the revised clouds to characterize the dis¬ 
tributions of the parameters. The parameters are listed 
in the supplement [3p| . It is important to note that this 
analysis is only possible because of the high precision 
data produced from the steady shocks. 

Figure [3] shows the compiled experimental Us — Up 
Hugoniot data, the four linear fits, and the phase regions 
determined from the MCO method. Following the lit¬ 
erature [25| and our DFT results, we propose the four 


regions be classified as follows: 1. The Bl solid from am¬ 
bient to 363 GPa; 2. The B2 solid from 363 to 462 GPa; 
3. The B2-liquid coexistence region between 462 and 
620 GPa; and 4. The liquid state above 620 GPa. How¬ 


ever, as our continuum level experiments do not provide 
microstructure information, we performed ab initio cal¬ 
culations of the Hugoniot and the phase diagram to bet¬ 
ter understand the high pressure states of MgO. 

The high precision requirements of this work necessi¬ 
tated refinements of previous ab initio methods 23-251]. 
We performed calculations utilizing DFT and QMC fo¬ 
cusing on the solid-solid phase transformation from Bl 
to B2 and the melting of MgO along the Hugoniot, pre¬ 
sumably from the B2 phase. Using DFT to calculate the 
Hugoniot requires prior knowledge of the phase state, so 
we first calculated the phase diagram. We used a three- 
part approach to determine the phase boundaries. In 
order to determine the melt boundary from both the Bl 
and B2 phases, we performed two-phase calculations of 
melting using VASP 5.2.1lE3, EH]i further details are 
presented in the supplemental material [30j. To deter¬ 
mine the solid-solid phase boundaries we decomposed the 
solid’s Helmholtz free energy into two pieces. 


F sol (V,T) = E(V) + F vib (V,T) 


( 1 ) 


The first piece is the density dependent energy of either 
the Bl or B2 phase. This is calculated via diffusion QMC 
usin g Q MCPACK [34J following methodology detailed in 
Ref. |35| with particular concern paid to the construction 
of pseudopotentials. The second piece of the free en¬ 
ergy is due to the finite temperature motion of the ions 
and electrons and is calculated in two parts. First the 
harmonic part of the free energy is calculated using the 
finite displacement method as implemented in the PHON 
code jHil] . The quasiharmonic approximation (QHA) is 
known to break down as temperatures increase and this 
is particularly true for MgO [37}. For this reason and be¬ 
cause the Hugoniot is expected to cross the phase bound¬ 
ary relatively close to the melt line, we have augmented 
our QHA calculations of free energy with thermodynamic 
integration. This is performed by using 

r T f 1 /dE 
T . T \dT 


AS = 


dT 


( 2 ) 


that allows the change in entropy along an isochore to be 
calculated directly in terms of the internal energy. The 
energy is calculated using DFT based quantum molecular 
dynamics (QMD) at points spaced by 250 K along sev¬ 
eral isochores in the region of the phase transition. Using 
entropy from the QHA calculation at low temperatures 
as a reference, we calculate the Gibbs free energy of both 
phases and determine the phase transition pressure di¬ 
rectly. This method also determines the range of validity 
for the QHA. We find the range to be smaller than pre¬ 
viously estimated US] with significant deviations in the 
free energy occurring by 5000 K and 400 GPa. The pos¬ 
itive effect of the anharmonic entropy was significantly 
larger in the Bl phase than in the B2 phase, moving the 
phase boundary to higher pressures at high temperature. 
Specific computational details are in the supplement E3 ■ 
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TABLE I. Phase boundaries along the principal Hugoniot. 


Method B1-B2 

B2-Coexist. Coexist.-Liquid 

(GPa) 

(GPa) 

(GPa) 

Expt. (MCO) 363±6 

462±20 

620±17 

DFT 330 

475 

600 


Figure 0] shows the DFT determined phase diagram 
and the calculated P-T states on the MgO Hugoniot. Us¬ 
ing the DFT phase diagram, we can infer the microstruc¬ 
ture during our plate impact experiments. However, the 
experimental Hugoniot data is an incomplete thermody¬ 
namic description due to the lack of temperature mea¬ 
surements. To assign the data to a given phase within 
the single phase regions, we use QMD calculations to pro¬ 
vide the missing temperature by calculating the Hugo¬ 
niot. For each of the three phases at given densities, 
we perform several QMD calculations at various tem¬ 
peratures. We then find the temperature for which the 
RH energy equation is satisfied to determine the shock 
state. Finally, the pressure and temperature of these 
shock states are compared to the phase boundaries to 
determine whether that state is thermodynamically sta¬ 
ble. This approach also provides a validation of the QMD 
by comparing the calculated Hugoniot to the experimen¬ 
tal results. The calculations and the experiments are in 
good agreement in p—P space (Fig. [2]) and in P—T space 
(Fig. 0|. The DFT calculated phase boundaries along the 
Hugoniot corroborate the MCO fitting method results for 
the experimental data suggesting the Hugoniot has four 
major regions: Bl, B2, coexistence, and liquid. Table Q] 
lists the phase boundaries along the principal Hugoniot 
from the MCO method and the DFT simulations. 

Both the experimental and DFT results show that 
a minimum shock pressure of 620 GPa is required to 
achieve complete melting in MgO that is initially at ambi¬ 
ent temperature. In the giant impact scenario, the proto- 
Earth has also been assumed to have an elevated surface 
temperature prior to the moon-forming event Q. We 
have performed additional DFT simulations to calculate 
the Hugoniot of MgO starting from an initial tempera¬ 
ture of 1900 K. From Tq = 1900K, a minimum shock 
pressure of 445 GPa is required to achieve complete melt 
in the MgO. Assuming planar normal impact, we can 
determine a minimum impact velocity required to melt 
MgO. Table ITT1 lists the required impact velocities for im- 
pactors of common planetary materials. In a real impact 
event, oblique impact (38] and shock attenuation through 
the mantle [39] affect the amount of energy transferred 
and imply even higher velocities are likely required for 
significant melting of MgO in the mantle. 

We have performed an extensive experimental and 
computational study of the high pressure - high temper¬ 
ature behavior of MgO. The Hugoniot is experimentally 
determined up to 1.2 TPa and the data suggests that the 


TABLE II. Impactor velocities for common planetary materi¬ 
als required to completely melt MgO with To at ambient and 
at 1900 K in a shock event (assuming planar normal impact). 


Initial MgO Temp. Impactor Impact Velocity 
[K] [300A'] [ km/s\ 

300 

MgO 

18.6 

300 

Dunite 

19.4 

300 

Iron 

15.3 

300 

Quartz 

20.1 

1900 

MgO 

16.0 

1900 

Dunite 

16.3 

1900 

Iron 

12.9 

1900 

Quartz 

17.7 



FIG. 4. P-T phase diagram of MgO with calculated Hugoniots 
starting at ambient and elevated initial temperature condi¬ 
tions. Experimental P-T data 02|, 00] and the low pressure 
B1-B2 melt line from Ref. FTil are included. 

Hugoniot crosses the B2 region with a large solid-liquid 
coexistence region lasting over 100 GPa before complete 
melt. The DFT simulation results corroborate the find¬ 
ings inferred from the experimental data and also show 
that at elevated T 0 the Hugoniot persists through a large 
coexistence region before melt. The results help place a 
lower bound on impact velocities for the moon-forming 
giant impact scenario. The data and phase diagram pro¬ 
vide a solid basis for the development of an MgO equation 
of state that can be used in planetary collision studies. 
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